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on 

ON ■ A variant of self-similar approximation theory is suggested, permitting an easy and accurate 

summation of divergent series consisting of only a few terms. The method is based on a power- 
law algebraic transformation, whose powers play the role of control functions governing the fastest 
convergence of the renormalized series. A striking relation between the theory of critical phenomena 
■ and optimal control theory is discovered: The critical indices are found to be directly related to 

limits of control functions at critical points. The method is applied to calculating the critical indices 
for several difficult problems. The results are in very good agreement with accurate numerical data. 
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The standard, and in many cases the sole, tool of theoretical physics is perturbation theory giving physical quantities 
in the form of power series. However, such perturbative series are usually divergent. If a sufficiently large number 
of terms of a divergent series are known, one can invoke resummation techniques. In theoretical physics, the most 
CO common resummation technique is Pade summation [1]. Unfortunately, Pade approximants and related techniques 
are often not able to produce sufficiently accurate results if only very few terms of a divergent series are available. 
Thus, one is frequently confronted with the annoying problem of extracting a sufficient amount of information from a 
divergent series even if only a few terms are known. An approach [2-6] , called the self-similar approximation theory, 
for solving this problem was developed, combining the renormalization-group techniques [2-4] with the methods of 
<—j . dynamical theory and optimal control theory [5,6]. One of the main ideas of this approach [2-6] is to introduce control 
functions governing the fastest convergence of approximation sequences [7] . It was shown recently [8] that an effective 
way of introducing control functions is by means of a power-law algebraic transform with powers playing the role 
of these control functions. Mathematical details of such an algebraic self-similar renormalization were expounded in 
7—i ! Ref.[8]. 

^ ■ In this paper we show that the approach can be extended to such a difficult problem as the calculation of critical 
indices for complex physical models, for which only a couple of perturbative terms can be obtained. A very interesting 
relation between control functions and critical indices is discovered. We illustrate the method by calculating the critical 
indices for several quite complicated systems. The results are in very good agreement with numerical data available. 
To make this paper self-consistent, we outline, first, the scheme of the method to be used. The complete mathemat- 
^ ical foundation can be found in earlier publications [2-8]. Since the most elaborate presentation of the approach [5,6,8] 
Q\ , relies on the language of dynamical theory and control theory, it would probably be helpful to give some additional 
■ references on dynamical theory [9] and optimal control theory [10] where a reader could find more details on the 
C]3 | mathematical language and techniques used in the approach. 

Assume that we are interested in a function tp(x) of a real variable x. Let perturbation theory give for this function 
perturbative approximations pk (x) with k = 0,1,2,... enumerating the approximation order. Define the algebraic 
transform Fk(x, c) = x c pk(x). This transform changes the powers of the series pk{x) changing by this the convergence 
properties of the latter. As a result, the approximation order effectively increases from k to fc+c. The inverse transform 
is pk{x) = x~ c Fk(x, c). Define the expansion function x = x(f, c) by the equation F (x, c) = /, where F is the first 
available approximation and / is a new variable. Substituting x(f, c) back to Fk, we get yk(f, c) = Fk(x(f, c), c). The 
. £h j transformation inverse to the latter reads Fk(x,c) — yk(Fo(x,c),c). 

Consider the family {?/&} as a dynamical system in discrete time. As far as the trajectory {yk(f, c)} of this 
dynamical system is, by construction, bijective to the approximation sequence {F)~(x, c)}, this system can be called 
the approximation cascade. Embed the discrete sequence {j/fe(/, c)} into a continuous sequence {y(r, /, c)} with 
t S [0, oo). Thence, the family {y(r, . . .) : re [0,oo)} composes a dynamical system with continuous time, whose 
trajectory passes through all points of the approximation cascade trajectory. Because of this, such a system can be 
called the approximation flow. 

The evolution equation for a flow can be presented in the functional form y(r + r', /, c) = y(r, y(r', /, c), c). We call 
this equation the self-similarity relation. Since the property of self-similarity is the central concept of our approach, 
let us remind where this notion is originated from. First of all, recall that this property is common for any autonomous 
dynamical system, reflecting the group property of motion [9]. Such relations appear in the renormalization-group 
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approach of quantum field theory [11]. Recently, the concept of self-similarity was widely employed in the context 
of fractals. In all the cases, self-similarity is the group property of a function conserving its form under the change 
of its variable. In particular cases, this variable can be time, as for dynamical systems, momentum, as for field 
theory, or space scale, as for fractals. In our case, such a variable is the approximation number, playing the role 
of time, and the motion occurs in the space of approximations, where self-similarity is a necessary condition for 
convergence [3,12]. The evolution equation for the approximation flow can be rewritten in the differential form 
and then integrated over time between k and some k*. The point k* is to be chosen so that to provide the best 
approximation F k+1 (x, c) = y(k*,F (x, c), c) for the minimal time k* — k. The cascade velocity Vk(y, c) in the vicinity 
of the time moment k may be presented by the Euler discretization of the flow velocity giving Vk(f, c) = V k (x(f, c), c), 
with Vk(x, c) — Fk+i{x, c) — F k (x, c). The integral form of the evolution equation is 

f F *+ l df 

JF k V k {f,C) 

where F k = F k (x, c) and F k+1 = Fj* +1 (x, c) . The approximation F k+1 must be reached during the minimal time. When 
no additional constraints are imposed, then the minimal time corresponds, evidently, to one step: k* = k + 1. Finding 
F k (x, c) from (1), and using the inverse transform, we come to the self-similar approximation p* k (x, c) — x~°F k (x, c). 

Now, by means of the substitution c — > c k , we have to introduce control functions c k which would govern the 
convergence of the sequence {p* k {x, c k )}. For the latter, following the standard procedure [2-6] and similar to the steps 
described above, we may construct an approximation cascade with the cascade velocity 

d 

V k *(x,c k ) = p* k+l {x,c k ) -p%(x,c k ) + (cfe+i - c k )—p* k (x,c k ), (2) 

which is the Euler discretization taking into account the finite-difference variation due to the varying control function. 
Convergence of an approximation sequence is, in the language of dynamical theory, the same as the existence of an 
attracting fixed point for the corresponding approximation cascade. If the cascade trajectory tends to a fixed point, 
this means that V k (x, c k ) — ► as k — ► oo. In practice, we cannot, of course, reach the limit k — > oo, and have to stop 
at a finite k. Then, the condition to be as close to a fixed point as possible is the minimum of the velocity: 

\V k *(x,c k (x)\ =mm\V k *(x,c)\. (3) 

c 

To make Eq.(3) simpler, we notice [4-6] that in the vicinity of a stable fixed point we have p k+1 ~ p k , c^+i w c k , and 
dp k /dc k sa 0. Consequently, the last term in the velocity (2) is of the second order of smallness, as compared to the 
difference p k+1 — p k - Therefore the practical way of using Eq.(3) is to start from the first-order term considering the 
minimal difference condition [7], that is, min c \pl +1 (x, c) — p k (x, c)|, and, only if this has several solutions, then one 
has to chose that one which minimizes the second-order term. Note that the fixed point equations as in Eq.(3) have 
been used in Ref.[13]. 

After the control functions are found from the minimal-velocity condition (3), we substitute them into pi and obtain 
the final expression f k (x) = p%(x, c k (x)) for the self-similar approximation of the sought function. To check whether 
the obtained sequence {f k {x)} converges, we have to analyse whether the corresponding mapping is contracting. The 
mapping related to the sequence {f k } is constructed in the standard way [2-8] and the contraction, or stability, is 
analysed by calculating the mapping multipliers. 

Now we shall illustrate by several examples how our method of algebraic self-similar renormalization works. We 
show that we can describe the critical behaviour well, starting from virial type expansions containing only second- 
order terms and derived for a region being far from the critical point. We would like to stress that the general approach 
formulated above naturally suits the complicated and unsolved problem of deriving the critical indices from very short 
virial-type series and that the control functions introduced by means of the algebraic transformation happen to be 
directly related to critical indices. Thus, it turns out surprisingly that critical indices can be considered as physical 
analogs of such mathematical objects as control functions. 

Begin with polymer coils. For these one usually considers the expansion factor a(z) as a function of the excluded- 
volume variable z = By/N(3/2irl 2 ) 3 / 2 , in which B is an effective binary cluster integral; N, the number of bonds 
of the length I each [14-16]. For the region of < z <C 1, corresponding to the poor-solvent or weak-coupling case, 
perturbation theory in powers of z gives [14-17] the expansion a 2 (z) ~ 1 + a\z + a 2 z 2 . On the other hand, in the 
good-solvent or strong-coupling case, when z»l, the polymer swelling is described by the law a 2 (z) ~ Az 2{ - 2lJ ~^\ 
in which v > 1/2 is a critical index. 

(i) 2D Polymer Coil. In the two-dimensional case, the coefficients in the virial-type perturbative expansion [15] 
are: a\ = 0.5 and a 2 = —0.12154525. Starting from the first-order and second-order perturbative expressions 
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Ph(z) for a 2 (z), following the prescribed way [8], as the self-similar approximations p£, when z — > oo, we obtain 
p1(z,c) ~ AfeZ~ c , where Ai = (z^r) arid ^2 = a i +c (zj^) • Comparing this with a 2 (z), as z — > 00, we get 



-01 

c = 2(1 — 2z/). Here c < is the limit, as z — > 00, of the control function c(z). Thus, we have the relation 

1 1 ,. \na(z) 11... 

1/=- + - hm - v ; = — — - hm c(z) (4) 

2 2 a-oo In z 2 4 z-oo w w 

between the critical index v and the limit of the control function c(z). From the fixed-point equation in the form 
of the minimal velocity condition (3), wc find c = —1, which gives v = 0.75. This is the exact value for the critical 
index [17]. Note that the Flory [18] mean-field formula v = 3/(2 + D), valid for D < 4, also reproduces this exact 
value. The reason why the mean-field approach can become exact can be understood if one notices [17] that the 2D 
polymer corresponding, within the framework of the 0(n)-model, to n — > belongs to the same line in the n — D 
plane as the Gaussian polymer with n = —2. Along this line, the mean-field gives exact critical indices [19]. 

(ii) 3D Polymer Coil. In the three-dimensional case, the coefficients of the perturbative expansion [14-16] are: 
a\ = 4/3 and a 2 = 287r/27 — 16/3. Now it is more convenient, for the stability of the procedure, to invert a 2 and to 
study a~ 2 (z) ~ 1 + b\z + b^z 2 , where b\ = —4/3 and 62 = 3.853. Repeating the same steps as for the 2-dimensional 
case, only with bk instead of a*;, we get c = 2(2f — 1). From the minimal velocity condition we find c = 0.396 , 
which yields v — 0.599. This is close to the Flory [18] and Edwards [20] mean-field result v = 0.6, as well as to the 
renormalization group calculations [21-24] giving v = 0.589. For the critical amplitude we find A = 1.62, which is 
also close to the renormalization group results [21-24] varying in the interval 1.53 < A < 1.75. 

The problem of the static conductivity of binary inhomogeneous materials is a well-known example of a rather 
complicated problem. Here one considers the resistivity p(x) and conductivity a(x) as functions of the concentration 
x of the superconducting or conducting bonds or sites in a disordered system [25,26]. The critical behaviour is 
characterized by the so-called superconductivity, s, and conductivity, t, exponents describing the power laws of 
resistivity p below the percolation threshold x c and of conductivity a above this threshold. The mean-field-type 
approximation for these disordered systems is the effective-medium approximation [25,27,28] which gives s = t = 1 
independently on space dimensionality. This is in sharp disagreement with the accurate numerical data [26] obtained 
by means of extensive transfer-matrix computer calculations. The position-space renormalization group [29-33] yields 
the results being also in a poor agreement with the numerical data [26] . 

(iii) 3D Disordered Superconductor. Below the percolation threshold, the effective-medium expression for the 
resistivity of bond-disordered systems [26-28] gives p(x) = 1 — 3x. From here we have for the self-similar renormalized 
resistivity p\{x, c) = c c /(c + 3x) c . This immediately gives us the relation 

s = lim ^P( x ) _ _ jj m c ( x ^ x < Xc; (5) 
x^x c m(x c — X) x^x c 

between the superconductivity exponent and the control function. From the equation p(x c ) = we find another 
relation, c = -3x c . With the threshold [26] x c = 0.2488, we obtain c = -0.746. Then from Eq.(5) we get s = 0.746, 
which is close to the numerical value s = 0.73 ± 0.01 from Rcf. [26]. 

(iv) 2D Disordered Conductor. Dealing with disordered conductors, it is convenient to exploit the fact that for 
them the effective-medium approximation becomes asymptotically exact in the low-defect-concentration limit, that 
is, in the pure-conductor limit x — > 1 [25-28]. Then, perturbation theory in powers of x — 1 gives [34] for the two- 
dimensional case a(x) ~ 1 — w(l — x) + 1.28588(1 — x) 2 . Following the same procedure of the algebraic self-similar 
renormalization, as is described above, we find the relation 

lner(x) 

t = hm — = — hm c(x), x > x c , (6) 

x^x c \n(x — x c j 



between the conductivity exponent and the control function. The latter is given by the minimal velocity criterion 
(3). The condition o-* k {x c ) = defines the threshold x c = 0.5905 which is close to the numerical value [26] being 
x c = 0.5927. From Eq.(6) we get t = 1.287, which is again in a good agreement with the numerical value [26] being 
t= 1.299 ±0.002. 

(v) 3D Disordered Conductor. For the three-dimensional cubic lattice, perturbation theory gives [25] the expansion 
a(x) ~ 1 — 2.52(1 — x) + 1.52(1 — x) 2 . From the threshold equation crl{x c ) =0we have x c = 0.255, which is close 
to the numerical result [26] for x c — 0.3116. The procedure of self-similar renormalization is the same as for the 
two dimensional case, yielding the same relation (6). For the conductivity index we find t = 1.88 which practically 
coincides with the numerical result [26] giving t = 1.9 ± 0.1. 
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To emphasize the generality of our method which works for perturbative series in powers of arbitrary expansion 
parameters, let us analyse one more nontrivial model. 

(vi) (2+l)D Ising Model. Consider the so-called (2 + l)-dimensional Ising model given by the hamiltonian H = 
J2i(^ ~ °f) — x J2<ij> <7 f cr | — a f > where i and j enumerate sites on the two-dimensional square lattice, (ij) 
denotes nearest-neighbour pairs, af are the Pauli matrices, x corresponds to the dimensionless inverse temperature 
in the Euclidian formulation, and h is the magnetic field variable. The high-temperature expansion [35] for the 
susceptibility at zero field yields x — 1 + 4a; + 13. 5x 2 , or for the inverse quantity ~ 1 — 4a; + 2.5x 2 . Following 
again the same procedure as above, we find the relation 

In v 

7 = - lim — = lim c(x), x < x c , (7) 

x^x c m(x c — X) x^x c 

between the critical index 7 and the control function c(x). For the critical point we get x c = 0.308 and for the index, 
7 = 1.232. This is in good agreement with the Pade summation of high-temperature series containing the terms up 
to sixteenth order [35], which gives x c = 0.3285 and 7 = 1.245 ± 0.005. 

In conclusion, we have advanced here a new method allowing an easy and accurate summation of divergent series 
containing only a few terms. A very interesting fact has been discovered that the critical indices are directly related 
to control functions governing the convergence of renormalized series and providing the stability of calculational 
procedure. The calculated values for the critical indices are in very good agreement with numerical data and exact 
results when the latter are available. The intriguing relation between the critical indices and control functions permits 
us to conjecture that the former play the same role in nature as the latter in theory: Nature choses such critical indices 
that provide the stability of physical phenomena in the vicinity of critical points. 
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